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ABSTRACT 

Context. Ogata et al. (2009; hereafter OKK) presented a theoretical determination of the ^Keiaa, yY^C, or triple-a, nuclear reaction 
rate. Their rate dilfers from the NACRE rate by many orders of magnitude at temperatures relevant for low mass stars. 
Aims. We explore the evolutionary implications of adopting the OKK triple-c reaction rate in low mass stars and compare the results 
with those obtained using the NACRE rate. 

Methods. The triple-a reaction rates are compared by following the evolution of stellar models at 1 and 1.5 with Z=0.0002 and 
Z=0.02. 

Results. Results show that the OKK rate has severe consequences for the late stages of stellar evolution in low mass stars. Most 
notable is the shortening-or disappearance-of the red giant phase. 

Conclusions. The OKK triple-or reaction rate is incompatible with observations of extended red giant branches and He burning stars 
in old stellar systems. 
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1. Introduction 

[ OKK introduced a new determination of the triple-ff nuclear re- 
. action rate based on a direct solution of the Schrodinger equa- 
' tion. The authors provided a tabulation of the triple-ff rate for 
temperatures from 10^ to 10'^ K and also the ratio of their rate to 
the NACRE (Angulo et al. 1999) rate. The OKK rate is a factor 
10^*= greater than the NACRE rate at 10^ K, lO*" times greater at 
10** K, and reaches equality at 25 x 10^ K and above. It is impor- 
tant to note that OKK normalized their rate to the NACRE result 
at 10^ K. For comparison, the previous update to the triple-ff 
reaction rate was from Fynbo et al. (2005). Between 10^ and 
10^ K, the rate determined by Fynbo et al. differed from the 
' NACRE rate by a factor of ~ 3 at most. The previous standard 
rate, that of Caughlan & Fowler (1988), was generally lower than 
the NACRE rate but never by more than a factor of ~ 10. 

Such a tremendous difference at low temperatures is certain 
to have dramatic implications for stellar evolution calculations 
and it was with this in mind that we undertook the following in- 
vestigation. The purpose of this research note is simple: to inves- 
tigate the evolutionary implications of the new triple-ff rate and 
compare results obtained from the new rate with results based 
on the NACRE rate. 

2. Stellar Evolution Model Comparison 

Stellar evolution calculations were carried out using two codes: 
DSEP (Dotter et al., 2007) and MESA (developed by Paxton and 
collaborators; http://mesa.sourceforge.net/). The input physics 
employed by DSEP and MESA are similar for the models pre- 
sented here. Both codes use a combination of OPAL (Iglesias 



& Rogers, 1996) and Ferguson et al. (2005) radiative opacities. 
For simplicity, the nuclear reaction networks in both codes were 
set to use the NACRE rates for all relevant H- and He-burning 
reactions, except that either the NACRE or OKK rate was used 
for the triple-ff procesfl For the equation of state, DSEP uses 
the FreeEOS code (http://freeeos.sourceforge.net) while MESA 
uses a combination of OPAL (Rogers & Nayfonov, 2002) and 
SCVH tables (Saumon et al. 1995) along with the HELM EOS 
(Timmes & Swesty, 2000) for temperatures and densities corre- 
sponding to He-burning and beyond. DSEP models use a mix- 
ing length parameter of umlt - 1-94 and MESA models use 
o'MLT = 2.14. The DSEP models include gravitational settling of 
He and metals while the MESA models do not. 

Stellar models were computed at 1 and 1.5 Mq with ini- 
tial X=0.70 and Z=0.02 using both codes. DSEP also produced 
models at the same masses with X=0.75, Z=0.0002. DSEP 
followed the evolution from the fully-convective pre-main se- 
quence phase until the onset of core He burning. MESA followed 
the evolution from the ZAMS through the core He flash and core 
He burning. The pre-main sequence evolution is not affected by 
the triple-ff process and is not shown in the following figures for 
clarity. 

Figure [T] shows the H-R diagram with DSEP models at 1 
and 1.5 M^: the left panel includes tracks at X=0.7,Z=0.02 and 
the right panel shows tracks at X=0.75,Z=0.0002. Both pan- 
els clearly demonstrate that the DSEP/NACRE models experi- 



' The original version of the OKK paper contained some typos in 
the tabulated rate. The correct rate, provided to D. VandenBerg by K. 
Ogata via private communication, was used in the present evolutionary 
sequences. 
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Fig. 1. DSEP models in the H-R diagram for X=0.7,Z=0.2 (left panel) and X=0.75,Z=0.0002 (right panel. Evolutionary tracks for 
1 and L5 are shown in both panels. Evolution from the ZAMS to the onset of the core He flash is shown in the 1 case. For 
the L5 Mq case with the OKK rate, the models transition to core He-burning without an explosive 'flash' event and a small portion 
of core He-burning evolution, which moves the model back to the blue, is shown. 




Log Tgjj, Age (Gyr) 

Fig. 2. The H-R (left panel) and age-luminosity diagrams from the 1 Mq, Z=0.02 MESA model. MESA is able to evolve the models 
through the core He flash and both tracks show core He-burning phases in the H-R diagram. The right panel demonstrates the 
prolonged core He-burning phase that results from adopting the OKK triple-a rate. 



ence an extended red giant phase in all cases. Conversely, the 
DSEP/OKK models show little or no red giant phase, with the 
Z=0.0002 models transitioning to core He-burning while still in 
the subgiant phase. The OKK rate also causes a very modest in- 
crease in surface temperature during the main sequence but the 
effect is small and difficult to see in the figure. 

Figure |2] shows, in the left panel, the H-R diagram of the 1 
Mq, Z=0.02 model evolved through the core He flash and core 
He-burning phase by the MESA code. The MESA/OKK model 



has a core He-burning phase that occurs at ~3-4 L^. about an 
order of magnitude fainter than the core He-burning phase ex- 
perienced by MESA/NACRE model. The right panel of Figure 
|2] shows the late-time evolution of the same models. Because 
the MESA/OKK model undergoes core He-burning at such low 
luminosity, it is able to continue for about 1 Gyr at more or 
less constant luminosity. The MESA/NACRE model experiences 
core He -burning for ~ 100 Myr at roughly ten times the luminos- 
ity of the MESA/OKK model. 
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Fig. 3. Central temperature-density diagram for 1 Mq models at Z=0.02. The left panel shows the DSEP models and labels the 
location of the ZAMS and core He flash for both models. The right panel shows MESA models, which also show the core He- 
burning phase of evolution. 



Figure [3] shows the central temperature-density diagrams for 
the 1 Mo, Z=0.02 models with DSEP on the left and MESA 
on the right. Models evolved with both codes trace out identi- 
cal paths in this plane until the DSEP models terminate at the 
onset of the core He flash. The MESA tracks show the evolution 
through the core He flash and core He-burning phases. The figure 
demonstrates that the OKK rate causes the core He flash to occur 
at a much lower temperature (-30 million K vs. -80 million K) 
and an order of magnitude lower density than the NACRE rate. 

3. Conclusions 

Comparisons of stellar evolution models computed with NACRE 
and OKK triple-a reaction rates were presented. The OKK rates 
cause models of low mass stars to have either a shortened red 
giant phase or to bypass the red giant phase altogether. The OKK 
rate also causes core He-burning to take place at a factor of ten 
lower luminosity-and for a factor ten longer duration-than the 
NACRE rate. 

Given the excellent agreement between stellar evolution 
models and observations of the RGB (Salaris, Cassisi, & Weiss, 
2002) when the models employ the NACRE (or relatively simi- 
lar) triple-a reaction rates, such as Caughlan & Fowler (1988) or 
Fynbo et al. (2005), it is difficult to envision how such a dramatic 
change as that proposed by OKK could lead to improved agree- 
ment with observations. Indeed, the triple-a reaction rate deter- 
mined by OKK is in stark disagreement with fundamental obser- 
vational evidence of stellar evolution, namely: (1) the existence 
of extended red giant branches in old stellar systems, (2) the lo- 
cation of core He-burning stars in the H-R or color-magnitude 
diagram, and (3) the lifetime ratio of the core He-burning phase 
to the red giant phase as manifested by the number ratio of hori- 
zontal branch stars to RGB stars, also known as the R parameter 
(e.g. Salaris et al. 2005). 
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